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ABSTRACT 



Propagation of the extremely short electromagnetic pulse in non-linear 
dielectric media without the slowly varying envelope approximation is dis- 
cussed. The models under consideration take into account both resonant 
and not-resonant excitations of non-linear medium, and polarisation states 
of electromagnetic wave. Steady state solutions of the relevant equations are 
presented for certain of these models. 
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1 Introduction 



In the context of rapid progress in the generation of femtosecond pulses of the 
electromagnetic radiation, it is of interest to consider the theoretical models 
of the pulses propagation in non-linear dispersive medium. Indeed, derivation 
of the evaluation equations for electromagnetic radiation is often based on 
the approximation of slowly varying complex envelopes of the optical pulses. 
This approximation is employed widely in theoretical studies devoted to non- 
linear coherent phenomena. However, it is of interest to find a method of 
describing extremely short pulse (ESP) evolution without using the slowly 
varying envelope approximation. 

A natural foundation for all the theories is the Maxwell equations supple- 
mented with equations governing the evolution of the polarisation or currents 
arising in a medium subject to electromagnetic radiation. Since it seldom 
occurs that exact analytical results can be obtained in such problems, var- 
ious approximations simplifying the consideration and allowing analytical 
expressions to be derived are often employed. 

The typical models of the non-linear medium are the ensembles of iV-level 
atoms or anharmonic oscillators. So, we can discuss the different degrees of 
approximation in the theory of interaction and propagation of optical pulses 
in framework of the resonant and non-resonant representation for medium 

2 Extremely short pulse propagation in res- 
onance medium 

Let the spectral width of ESP be far smaller than the resonance frequency. 
This resonant condition allows to omit all non-resonant levels and to consider 
only two levels. Thus, we obtain the approximation of the two-level atoms 
|2|, which is very popular in the resonant and coherent optics. The wave 
equation in this case is complemented by the Bloch equations for two-level 
atom variables. 

2.1 Propagation of the scalar ESP 

In general, the theory of the interaction of radiation with an ensemble of 
two-level atoms is based on the Bloch equations for atoms and the Maxwell 
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equations for the classical electromagnetic field. In an isotropic dielectric 
the set of Maxwell equations reduced to one equation for the electric field 

— # — * 

E = El. For a plane wave with constant polarisation vector one can obtain 
the following system of total Maxwell-Bloch (MB) equations 



d 2 E 1 d 2 E A7m A d I d 2 r l 



dz 2 c 2 dt 2 c 2 \ dt 2 



(1) 



dr x dr 2 2d dr 3 2d 

^t=-^ ^t =UJari + T Er3 > ^t=~T Er2 > (2) 

where d is the projection of a matrix element of the dipole operator on the 
direction of I, n A is the concentration of resonant atoms. It should be noted 
that the components of Bloch vector r\, r 2 , and r 3 depend on the atomic 
resonance frequency io a . Hereafter the angular brackets represent summation 
over all the atoms characterised by the frequency uo a . 

Let us consider the solution of equations ([[]) and (0) describing a solitary 
steady-state wave. Boundary conditions at \t\ — > oo must look like the 
following 

E = dE/dt = 0, r 2 = n = 0, r 3 = -1. 

To obtain the expressions describing the propagation of stationary ESP 
one should suppose that the components of the Bloch vector and the electric 
field E depend on only one variable t±z/V or r] = u a (t±z/V) in dimension- 
less form. That means that a stationary (steady state) wave propagates only 
in one (positive or negative) direction. Under this assumption the system 
([]]), (fj) transforms into the system of ordinary differential equations which 
can be solved by usual methods. Solution of the these ordinary differential 
equations leads to the expression for electric field strength 

E(z,t) = E sech[dE (t±z/V)/h), (3) 

where pulse amplitude is related with pulse duration t p as E = ht~ x d~ x . 
This expression for stationary pulse correlates with the one found in 0, 
Expression for the stationary ESP velocity can be obtained from definitions 
of pulse duration and amplitude. As in self-induced theory by McCall-Hahn 
|2|. |5J, these values are related to each other, so the pulse is able to invert a 
two-level system and then return it to the initial state while the pulse lasts. 
After some algebra one may find 
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(4) 
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that corresponds to what was found in ||. 



2.2 Propagation of the ESP of polarised radiation 

In this case generalised Maxwell-Bloch equations can be written as follows 

d 2 E (+l) y Q2 E (+l) AnnA Q2 



dz 2 c 2 dt 2 c 2 dt 2 

d 2 E (-l) l Q2 E (-l) _ 4nnA Q2 

dz 2 c 2 dt 2 c 2 dt 2 



(d 13 p 31 + d 31 p l3 ) , (5) 

(<W 3 2 + d 32 p 23 ) , (6) 



= + dl3( ^3 " Pu)E {+l) ~ d23P 12 E { - 1 \ 

iH ^W = ~^ aP23 + rf23(/?33 ~ P^ E( ~ 1} - di3p2iE {+l \ 

ih 9 -^ = d 13 p 32 E^ - d 32Pl3 E^\ (7) 

d 

ih fo(Pn - P33) = 2 ( rf 13P 3 l - rf 3lPl 3 )^ (+1) - ( rf 2 3 P 3 2 - d32p23) Ei ~ 1] i 



d 

ih Ql(p22 - P33) = (<W 3 1 - d31Pi 3 )E {+1} - 2(d 2 3p 3 2 ~ d 32 p 23 ) E^ 1 ' , 

where is the spherical q component of the vector of electrical field 
strength of a light wave, q = ±1, dki are the matrix elements of the dipole 
moment operator of the atomic transition. Hereafter d\ 3 = c?23 = d 3l = ^32- 
By supposing that the matrix elements of the density matrix p kl , and the 
electric field E depend on only one variable t±z/V we can find the expression 
for the strength of the electrical field of the ESP in the form 

E {±1) (z,t) = e {±) E sec h[dE (t±z/V)/h] (8) 

where are the components of a unit vector of electromagnetic fields. This 
is a simple generalisation of the results obtained by Bullough and Ahmad 
0] to the case of vector (polarised) ESP and specific model of a resonant 
medium. As in Ref. JIJ, the pulse duration tp is expressed in terms of the 
pulse peak amplitude t p = h^dE^)' 1 . The propagation velocity of the ESP 
is: 

1 1 ( ^ 8nn A d 2 hu a \ ^ 



V 2 c 2 \ 2(dE ) 2 + (hu a ) 
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Besides the solitary wave solutions, equations ([I]), (0) and (|3|)-(0) have 
the solutions in the form of cnoidal wave. 

The above analysis of the ESP propagation shows that the nature of such 
stationary ESP depends on the state of the medium. The scalar solution |4| 
can be trivially generalised to cover a vector video pulse case under condition 
n w = n 2o = 1- So we obtain a circular polarised pulse with duration of 
a half of a reciprocal atomic transition frequency. The propagation velocity 
of this pulse does not depend on the conditions of its polarisation and the 
value of this velocity is in agreement with what was found in [4j]. 

The new solution of the generalised complete Maxwell-Bloch equations 
arises when initially the resonance medium has an asymmetrical distribution 
of level populations (i.e., n w and n 2 o) of the excited states that have different 
projections of the angular momentum: n\n ^ n 2 o 0. One of the spherical 
components of the electric field strength vector behaves as in the scalar case, 
i.e., it is a unipolar spike of the electrical field. Another component is a 
variable-sign solitary wave. All stationary solutions like these form a three- 
parametric family, where each member is defined by the propagation speed 
V and n w , n 20 determining the polarisation state of the pulse. 

3 Extremely short pulse propagation under 
quasi-resonance condition 

Let assume that the amplitude of electrical field of electromagnetic pulse is 
so small that the instant Rabi frequency uir proves to be far less then the 
frequency of the resonance transition, i.e., 

e = — = 1 < 1, (10) 

where ujr corresponds to the wave of constant electrical field amplitude E = 
max \E^\. In this case the solution of Bloch equations ( |7|) can be presented 
formally || by power series in e. Limiting ourselves to the terms of this series 
up to e 3 , for polarisation in equations (|) and (||), p( +1 ) = c/ 13 p 31 + d 31 p 13 , 
p(~V = d 2 3P 32 + ^32P23> we obtain expressions [|J : 
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The substitution of (O) into equations (|5[) and (|G[) results in the non-linear 
wave equation, which gives approximate description of the extremely-short 
pulse of electromagnetic wave 

d 2 E 1 d 2 E d 2 r/ - 

-b^r-a^{(E-E)E} = (12) 



where 



dz 2 V 2 dt 2 dt A dt 2 

1 16iinAO'\di 3 \' 2 \ / 8^^(71^13 1 



\ c 2 (/iu; a ) 3 / ' \ c 2 ^ a 3 

The velocity of pulse propagation V alters because of the dispersion, which 
is due to the resonant medium: 

V~ 2 = c' 2 [1 - (8im A a\d 13 \ 2 /hio a )] . 

For zero boundary conditions as < -> ±oo the steady state solution of 
equation (|12"D can be written in the form: 

E(z,t) = e£ sech[(t ±az/V -t )/t s ] (13) 

where 



t s = £ ~Vl 6 l/ 2 R a = yjl + 2aE 2 V 2 \a\. 

Thus, the steady state solutions of the equation flT2"D describe the ex- 
tremely short pulse propagation with constant polarisation state (i.e., e = 
const ) of this electromagnetic wave. 

Spectral half- width of the ESP (^) is tg . Consequently the resonance 
condition will be applicable if tg 1 to. From the definition of the parameters 
\a\ and |6| we can estimate their ratio as ~ 2\d 13 \ 2 uj 2 h~ 2 (where we 

neglected inhomogeneous broadening and all the frequencies ui a were replaced 
by u). Thus, the resonance condition complies with the condition ( |T0D uor 
Co. 

Besides the solitary steady state wave ( |T3D both dark solitary waves and 
a periodic stationary waves can be found. However they contradict with 
conditions of applicability of the initial equations for atomic variables, in 
where any relaxation processes are not taking into account. The considering 
pulse is short one in this sense. 
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4 Extremely short pulse propagation in non- 
resonance medium 



As it was pointed in || when the ESP field approaches the atomic field (~ 10 8 
- 10 9 V/cm), ESP formation is affected mainly by the atomic ionisation 
potential, which limits pulse duration and amplitude. For such fields the 
resonant models are invalid. However, one can evaluate limitations on ESP 
within a classical model of an atom with a strong non-linear potential, limited 
at infinity, to allow for ionisation. Then the (generalised) Bloch equations 
are replaced by classical equations for the electron (and ion) motion. 



4.1 Scalar Duffing model 

Let the medium be a non-linear and non-resonant one. A similar approach 
was exploited [HJ to develop the self-induced transparency in ionic crystals 
in the framework of the Dufling's type model. The scalar wave equation is 

d 2 E ld 2 E_4nd 2 P 

We complete this wave equation by the non-linear oscillator equation for 
the medium. If X is the displacement of an electron from its equilibrium 
position, the motion equation (which neglects friction) can be written as 

^ + u 2 Q X + k 3 X 3 = Z-E (15) 

where uq is an eigenfrequency of the oscillator, while K3 is anharmonicity 
coefficient. The term on the right-hand side of equation (|15|) represents the 
force exerted on the electron by the electromagnetic field, where £ = (e+2)/3 
is the Lorentz factor, e is the electron's electric charge, and m is its mass. 
However, it is possible to absorb £ into an effective mass m e ff = m/£. 
Hereafter, we will use m as a symbol for this effective mass. Finally, the 
dynamical variable X is related to the medium's polarisation, P = riAeX, 
where tia is the density of the oscillators (atoms). 

The steady state solutions of the system of equations ( |14D and (|T5|) can be 
found if we suppose that the displacement X and electric strength E depend 
on one variable rj = u>o(t ± z/V). Here velocity V is a parameter of the 
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solution. Thus, the expression for electric field of ESP can be obtained by 
solving the ordinary non-linear differential equations. That leads to 

E(t ± z/V) = E m a(a - 1) 1/2 sec h [(a - l) 1/2 u (t ± z/V - to)] (16) 

where E m = {mu'l/ e) 2 / k^,, a = (uj p /ujo) 2 V 2 /(c 2 — V 2 ), andu; p = (AniiAe 2 /m) 1 / 2 
is the plasma frequency: 



4.2 Vector Duffing model 

Let us consider the ESP propagation in a non-linear medium, which is again 
presented by an ensemble of the molecules. However, now let us suppose 
that the inner degrees of freedom are described by the potential: 

U(X, Y) = ^X 2 + ^oo 2 Y 2 + X -k 2 X 2 Y 2 + ^ 4 (X 4 + y 4 ), (17) 

where k 2 and K4 are anharmonicity coefficients. If the coupling parameter 
k 2 = 0, then this potential describes the scalar Duffing model of anharmonic 
oscillator, which can oscillate independently in two orthogonal directions. 
Expression ( |17D is the simplest generalisation of this model. Polarisation of 
the molecule is defined by the expression p = eXei + eYe 2 , and the total 
polarisation P is the product of the density ha and the polarisation of one 
molecule. 

The propagation of the ESP will be considered in the unidirectional wave 
approximation, so the wave equations for electric field corresponding to the 
different polarisation components of the ESP can be written as follows 

d 2 E 1 1 d 2 E 1 _ Aim A e d 2 X d 2 E 2 1 d 2 E 2 _ Aim A e d 2 Y 

~dz T ~ ~d*~dt^ ~ c 2 dt 2 ' Ih? ~ ^~d^ ~ c 2 dt 2 ' ^ 18 ^ 

The equations of motion for anharmonic oscillator follow from the classical 
Newton equations 

^ + u\X + k 2 XY 2 + k 4 X 3 = -E 1: (19) 

+ lu 2 2 Y + k 2 YX 2 + k a Y z = -E 2 . (20) 
ot 2 m 

It is possible to find a steady state solution of this system of equations. 

Let us suppose the fields depend only on r = t± z/V . Integrating equations 
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and taking account of the boundary conditions, which correspond to the 
vanishing electric field of the USP and polarisation of molecules at infinity, 
we find that E\ = A^ua^OX and E 2 = AimAeOY. Substitution of this result 
into (|T9|) and ( p0|) leads to the following equations 



d 2 X 
'JP 2 

d 2 Y 
dr 2 



+ (k 2 Y 2 + k±X 2 )X = (9lu 2 - ui)X, (21) 



+ (k 2 X 2 + k 4 Y 2 )Y = {Out - oo 2 )X, (22) 



where 9 = V 2 (c 2 — V 2 )~ l These equations appear in many works, so we can 
easily find their solutions. Let suppose that the eigenfrequencies uii and uj 2 
are equal. The partial solution of system ([H]), (^) in this case is 

Ei(t) = E 2 {t) = E m a{a - l) 1/2 sech[(a - l) 1 ' 2 ^ - r )], 

where E m = {muj\/e)\j2/ (k 2 + k 4 ), t is the integration constant. This so- 
lution similar to that above describes the linear polarised video pulse prop- 
agating without distortion in the considered non-linear medium. 

The more interesting solution exists if the frequencies ui\ and uj 2 are differ- 
ent, but factors of anharmonicity coincide k 2 = k 4 . In this case the solution 
of the system of equations fl21 ) and (p2|) leads to following expressions JTT 



Ei(t) 



E 2 (t) 



where 



87ren^6 , y / 2fi^ T /i 1 exp(#i){l + exp(26* 2 + ^12)} 
1 + exp(20i) + exp(2# 2 ) + exp(20i + 29 2 + ^ 12 ) : 

8nenA0 \/%k 2 T /J> 2 exp(# 2 ){l + exp(29 1 + /i 12 )} 
1 + exp(20i) + exp(2# 2 ) + exp(20i + 26 2 + n l2 ) ' 

exp n 12 = {ii x - fi 2 ) {n x + fi 2 ) - 1 . 

2 , ,2 



(23) 



In these expressions the following parameters /i 12 = (9u 2 — u{ 2 ) and #i j2 = 
/i 12 (r — ri j2 ) are used, where ti j2 are the integration constants. Solutions 
(|23|) describe the steady state propagation of the ESP where one of these 
components corresponds to a unipolar spike of the electric field, sub-cycle 
pulse. The other component is a bipolar impulse or monocycle pulse. 
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4.3 Cubic-quadratic Duffing model 

An objective of the present section is to discuss the propagation and interac- 
tions of linearly polarised USPs in a non-linear dispersive medium modelled 
by an anharmonic oscillator combining quadratic and cubic nonlinearities. As 
is well known, in the case when the oscillator is quasi-harmonic, the quadratic 
and cubic non-linear terms produce effects of the same order of magnitude 



12| , that is why it is natural to consider a mixed model of this type. It is 
also relevant to mention that modelling dynamics of the broad (rather than 
ultra-short) optical solitons in a medium with competing quadratic and cu- 
bic nonlinearities has recently attracted considerable attention, see, e.g., the 
works |13|] and references therein. 



The one- dimensional propagation of the electromagnetic waves in a non- 



linear medium is governed by the wave equation (0). We adopt a simple 
anharmonic oscillator model for the medium, which is a frequently used ap- 
proximation [rj] (see also [T3, |HJ). If X is the displacement of an electron 



from its equilibrium position, the motion equation (which neglects friction) 
can be written as 

* + - K2 X 2 + k 3 X 3 = —E, (24) 
at 2 m 

where k 2 and are anharmonicity coefficients. It seems plausible that the 
system of equations ([14]), ( [24]) is not an integrable one. Nevertheless, some 
exact analytical solutions, describing the propagation of ESPs, can be found. 
To this end, one assumes that E and X depend on a single variable, rj = 
u)o(t ± z/V) with some velocity V. It is convenient to express the velocity V 
of a steadily moving pulse by introducing the constant a 

where ui p is the plasma frequency. The first equation of the system, i.e., 
([14]) , can be integrated. It allow us to obtain a family of exact solutions 
parameterised by the continuous positive parameter (a — 1) and discrete one 
a = ±1, 

E'^V, ») = 3£ - a( °" 1) , . (26) 

y/l + 9{a - l)/icosh f y/(a - 1)77 J - 1 



10 



&~\ m a) = 3E m a (a-l) ^ 

a/1 + 9 (a - Tj/1 cosh (•>/(« - !)??) + 1 
where £7 m = (mo^/el/^!), A* = (^3^0/2^2); an d the superscript standing 



for cr. The pulses represented by the solutions (|33|) and fl34|) have different 
polarities. In the limit K3 — > 0, corresponding to the model with a purely 
quadratic nonlinearity, equation (E6f) goes over into a singular solution, 



2sinh 2 (y(a- 1)77/2 



while expression ( p7| ) yields a non-singular pulse in the same limit, which was 
already found in Ref. [|17| , 



2 cosh 2 [^/{a- l)?7/2 

Thus, in compliance with the results of Ref. [|l7j, only one family of non- 
singular pulses exists in the case of the purely quadratic nonlinearity. We 
stress that the velocities of two pulses, which have the opposite polarities but 
the same value of a are equal, since they are determined by equation Q25| ) 
and depend only on a. 

Although the pulse solutions have been obtained in the exact form, their 
stability and interactions should be studied by means of numerical simu- 
lations. This investigation was performed under unidirectional approxima- 
tion. Simulations of interactions between have shown that their collisions 
are quasi-elastic, irrespective of the polarities of the colliding pulses: after 
passing through each other, the pulses retrieve the same shapes and veloc- 
ities as they had before the collision, the only result being a shift of their 
centers. However, the character of the interaction between pulses with equal 
and opposite polarities becomes different with the decrease of their relative 
velocity. In the case of equal polarities, the pulses with a small relative 
velocity demonstrate strong mutual repulsion: the distance between them 
attains a minimum, and then they start to separate again, so that they never 
completely overlap. Strong energy exchange between the pulses takes place 
around the point where they attain the minimum separation. The energy 
exchange gives rise to mutual interconversion of the two pulses, so that after 
the collision they, effectively, swap their positions. 
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4.4 Vector model of an quadratic non-linear media 



The propagation of an ESP of polarised radiation in a medium with a 
quadratic nonlinearity can be considered by using a model of an ensemble 
of two-component anharmonic oscillators . This model described by the 
wave equations fll8|) supplemented with the following equations 



* + U 2 X + Klll X 2 + 2k U2 XY + k 122 Y 2 = —E 1} (28) 
ot mi 

d 2 Y 



... + uoiY + Kn 2 X 2 + 2k U2 XY + k 222 Y 2 = —E 2} (29) 
Ot m J < 2 

where the anharmonicity coefficients k^i and two effective masses mi and m 2 
are introduced. Let the electric field strengths and displacements X and Y 
depend only on the variable r = t ± z/V. Similar to the scalar case, one can 
integrate the wave equations ( p~8|) and equations (p8[), (p9|) can be reduced 
into the following set of non-linear equations 

>b A + {uj\ - 9uj 2 p )X + KmX 2 + 2ku 2 XY + k 122 Y 2 = 0, (30) 



dr 2 
d 2 Y 



+ {u 2 2 - Bul)Y + ku 2 X 2 + 2k U2 XY + k 222 Y 2 = 0. (31) 



If both the eigenfrequencies of oscillators and the anharmonicity coefficients 
are equal to each other, then the problem considered can by reduced to the 
scalar one. Unfortunately, the general solution of the equations fl3Tjp, Q3"T| ) is 
unknown. 



4.5 Coupled between vibrations and light waves in Ra- 
man media 

To description of the light propagation in Raman medium in [^] the oscillator 
model has been proposed. The displacement Q, vibration density p, and 
electric field of the light wave E are governed by the following equations 

8 2 E 1 8 2 E 4tt 8 2 P 



dz 2 c 2 dt 2 c 2 dt 2 ' 

d 2 Q , 1 dQ 2 at , 

+ u oQ = -7TT7 E P, ( 32 ) 



dt 2 T 2 dt u ^ 2M 
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dp i «i 2 ag 

where u is an eigenfrequency of the oscillator, M is a mass of ion. The non- 
linear polarisation is defined by expression P = uaO-iQE, where parameter 
a,\ is the coefficient in expansion of the atomic (or molecular) polarizability 
a(Q) ~ ao + a\Q. The relaxation processes were taken into account by 
introduction relaxation times 7\ and T 2 . The constraint p = const = p = 
— 1 leads this model to the Placzek one describing the stimulated Raman 
scattering |21J. The typical magnitude of the parameters are ujq = 10 13 s _1 , 



M = 10~ 22 g, on = 10~ 15 cm 2 , n A \p \ = I0 22 cm~' 6 . 

The steady state solutions of the Blomberge-Shen model (^) are not 
known. It was shown that under unidirectional approximation (when wave 



propagates only in one of the possible directions) [21] this model can be used 



to description of the self-induced transparency under two-photon resonance. 



4.6 Isotropic solid dielectric model 



As pointed out in p|j , in the study of the ESP propagation in a solid dielec- 
tric medium it is essential that the material equations adequately account for 
dispersion of both the linear and the non-linear susceptibility of the medium 
over practically its entire transparency range. Taking into account the con- 
tributions of both the electron and phonon subsystems of the solid dielectric 
one can describe normal and anomalous group dispersion. Thus, the simplest 
generation of the Lorentz model was proposed 

^ + (u\ + 7 E 2 )P = n A a E + PRE, ^ + u\R = X PE, (33) 

where P is electronic polarisation, the extra variable R reflects vibration 
motion of the atoms. Parameters ao, /3, 7, x are the characteristic ones to 
phenomenological description of the dielectric medium. The more complete 
model was proposed in p3|. The generalised Lorentz model represented there 
by the set of equations 

fpp 9 BP 

~W + T~at + uj2ePe = aeE + p{Re + Rv)E ' 

d 2 P t 2 BP 2 



+ Tz-xr+vfPi = aiE 



Bt 2 Tt Bt 
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where P e and Pi are the contributions to the polarisation of the electronic and 
ionic components. The oscillators described by the variables R e and R v give 
rise to non-linear parametric coupling between the ESP field and the medium. 
The dynamical parameter R e is responsible for the electronic nonlinearity, 
and the dynamical parameter R v reflects electron-phonon nonlinearity. The 
phenomenological parameters a e>i , u e>i , T e>i characterise the dispersion of 
the electronic and ionic linear polarisation responses. The coefficients (3, 
j ev , u> e i jV , T e i jV characterise the dispersion of the electronic and electron- 
vibrational non-linear responses, respectively. 

A sequence of the approximations (it has been detailed in p3| ) allows to 
solve the equations (^) that leads to non-linear wave equation 

T T 

^ - + b J Edr' + gE^ = A ± J Edr', (35) 

— oo — oo 

where t = t — z/V is retarded time, a, b, g, and $ are positive parameters 
determined in ||23|1 . is the transverse Laplacian. The Ref. [p4J] represents 
the vector generation of the non-linear wave equation 

dE <9 3 £ ./"-., , - dE - , dE 



a- 



dz dr 3 



EdT' + g(E-E)— + hEx(Ex—)=0, (36) 



describing propagation of the ESP of arbitrary polarisation. The numer- 
ical simulations shown that the polarisation self-action of ESP consists in 
a variation of the orientation of the vector of the electric field strength in 
the direction to the propagation axis proportional to the field squared and 
the velocity of rotation of its strength vector. There are no exact analytical 
results, unfortunately. 

4.7 Gas of quantum particles 

A self-consistent model for a non-linear interaction of an intense extremely 
short light pulse with a gas of atoms in framework of the quantum mechan- 



ics has been discussed in [25|. The model includes the scalar wave equation 
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(TJ) and one- dimensional Schrodinger equation for an electron, which is sub- 
jected to the intratomic potential U(x) and the optical radiation in the dipole 
approximation 

dip h 2 d 2 ip , 
lh M + 2^d^- U{x ^ = exE ^ 
The polarisation P of the medium, which consists of model one-electron 
atoms with a number density ha, is 

+00 

P(t) = —eriA I x\ip\ 2 dx. 



By using the Schrodinger equation we can rewrite the wave equation as 

+00 

d 2 E 1 d 2 E 47re 2 n, 4vren A f . , . .. 2 dU 



dz 2 c 2 dt 2 c 2 m 



S+^/lWlffc (37) 
c z m J ox 



The first term describes the linear response of the free electrons in the ab- 
sence of the intratomic confining potential. The second term reflects the 
contribution of the coupled electrons. In |25] it was shown by numerical 



method that it is posible to generate a burst of harmonics with a relative 
power (in the region of 1%) and a duration amounting toonly a few periods 
of the exciting optical field. 

4.8 Non-truncated Maxwell-Bloch model 

In order to investigate the self-induced transparency phenomena for ex- 
tremely short electromagnetic pulses in ^Zfi\ the non-truncated Maxwell-Bloch 



equations have been proposed. These equations take the form 
d 2 A Id 2 A 4vr . A . ... 2 . 

d ^ t -i{u + 1 \A\ 2 )3 = i^fpA, (38) 

where A is a vector-potential of the electromagnetic wave, j is atomic transi- 
tion current density, p is a population inverse density, mi2 is matrix element 
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of resonant transition current. The Kerr and Stark effects were taken into 
account by introduction of the terms containing parameters X3 an d 7 S , re- 
spectively. 

The results of the numerical simulations demonstrate the existence of 
the self-similar solutions of the system of equations fl38f ). These solutions 
should be called by the solitons because they are stable in respect with pulse 
collisions. The main distinctive feature of the solitons of the non-truncated 
M axwell-Bloch equations in that they do not break-up into the pulse sequence 
when the area of the incident pulse is higher then 2n. 

5 Conclusion 

We have considered some case of model media that allow the explicit temporal 
dependence of the electric field strength of the extremely short pulses to be 
determined by analytically. We did not make any assumptions concerning 
the harmonic carrier of the wave or the variation rate of the field in the pulse. 
As usually, only steady state solutions describing the propagation of the ESP 
can be found. 

There are some investigations in the considered area that were left aside. 
However, it seems useful to mention works |27|, |2S|, [2D|] where the propagation 
of a powerful electromagnetic pulses in a strongly non-linear medium without 
dispersion has been considered. The exactly solvable models are developed to 
describe the interaction of extremely short (few-cycle) transients with certain 
classes of insulators and conductors. Transient-excited fields are described 
analytically based on new, exact, non-periodic and non-stationary solutions 
to Maxwell's equations, obtained directly in time domain using a no Fourier- 
expansion, no time-space separation method. Such non-separable solutions 
form the mathematical basis of the non-periodic electromagnetic waves and 
make the standard harmonic wave concepts of frequency, phase, refraction 
index, and phase velocity irrelevant to the time description of non-periodic 
waves. Ref. J^J] reviews these investigations. 
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